The close T Tauri binary system V4046 Sgr: 
Rotationally modulated X-ray emission from accretion 

shocks 



C. Argiroffi 1 ' 2 , A. Maggio 2 , T. Montmerle 3 , D. P. Huenemoerder 4 , 
E. Alecian 5 , M. Audard 6 ' 7 , J. Bouvier 8 , F. Damiani 2 , J.-F. Donati 9 , 
S. G. Gregory 10 , M. Giidel 11 , G. A. J. Hussain 12 , J. H. Kastner 13 , and 

G. G. Sacco 13 



ABSTRACT 

We report initial results from a quasi-simultaneous X-ray/optical observing 
campaign targeting V4046 Sgr, a close, synchronous-rotating classical T Tauri 
star (CTTS) binary in which both components are actively accreting. V4046 Sgr 
is a strong X-ray source, with the X-rays mainly arising from high-density 
(n e ~ 10 n ~ 12 cm~ 3 ) plasma at temperatures of 3 — 4MK. Our multiwavelength 
campaign aims to simultaneously constrain the properties of this X-ray emitting 
plasma, the large scale magnetic field, and the accretion geometry. In this paper, 
we present key results obtained via time-resolved X-ray grating spectra, gathered 
in a 360 ks XMM-Newton observation that covered 2.2 system rotations. We find 
that the emission lines produced by this high-density plasma display periodic 
flux variations with a measured period, 1.22 ± 0.01 d, that is precisely half that 
of the binary star system (2.42 d). The observed rotational modulation can be 
explained assuming that the high-density plasma occupies small portions of the 
stellar surfaces, corotating with the stars, and that the high-density plasma is 
not azimuthally symmetrically distributed with respect to the rotational axis of 
each star. These results strongly support models in which high-density, X-ray- 
emitting CTTS plasma is material heated in accretion shocks, located at the base 
of accretion flows tied to the system by magnetic field lines. 

Subject headings: Accretion, accretion disks — Stars: individual: V4046 Sgr — Stars: 
magnetic field — Stars: pre-main sequence — Stars: variables: T Tauri, Herbig Ae/Be 
— X-rays: stars 
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1. Introduction 

In the context of star formation and evo- 
lution, understanding the physics of young 
low-mass stars is essential. Such stars 
possess strong magnetic fields that reg- 
ulate the transfer of mass and angular 
momentum to and from the circumstellar 
disk, via accretion and outflow phenom- 
ena. Young low-mass stars are also intense 
sources of high-energy emission (UV and 
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X-rays) that ionizes, heats, and photoevap- 
orates material in the circumstellar disk, 
thus affecting its physical and chemical 
evolution and, eventually, the disk lifetime 
(lErcolano et al.l 120081 : iGorti fc Hollenbach 



2009|). 

Low-mass pre-main sequence stars are 
classified as classical T Tauri stars (CTTS) 
when they still accrete mass from the 
circumstellar disk. They become weak- 
line T Tauri stars (WTTS) when the ac- 
cretion process ends. Both CTTS and 
WTTS are bright in X-rays due to the 
presence of hot coronal plasmas, heated 



netic fields ( 


Feieelson & Montmerlel 


L999|; 


Favata & Micelall2003l 


Preibisch et al. 


2005 


Giidel & Naze! 


2009). 


It was suggested 



that in CTTS also the accretion process, 
beside the coronal magnetic activity, can 
provi de a further X-ray emission mecha- 
nism ( )Ulrichlll976t lGullbringlll994j ; lLamzin 
19991 ). Magnetospheric accretion models 



predict that in CTTS mass transfer from 
the inner disk onto the star occurs via ac- 



tubes 



(e.g. 


Koniel 


1991; 


Hartmann et al. 


Bouvier et al. 


2007), where material 



moves in a almost free fall with typical 
velocities of ~ 300 — 500 km s -1 . The im- 
pact with the stellar atmosphere, usually 
involving small fractions of the stellar sur- 
face, generates shock fronts that heat the 
infalling material up to temperatures of a 
few MK, and therefore should yield sig- 
nificant emission in the soft X-ray band 
(0.1 — 1 keV). Numerical modeling predicts 
high Lx (~ lO^ergs" 1 ) even for low ac- 
cretion rates (10~ 10 M© yr _1 ), indicating 
that X-ray emission related to the accre- 
tion process can rival or exceed coronal 



2 



emiss ion ( jGiinther et al.ll2007t ISacco et al. 



20081 ) at least in principle. 

Strong evidence of accretion-driven X- 
rays from CTTS has been provided by 
the observed high densities of the X- 
rav emitting plasma at T ~ 2 — 4MK 



~ 10 12 - 10 13 cm" 3 , 


Kastner et al. 




Schmitt et al. 


20051 


Giinther et al. 


Argiroffi et al. 


2007 


Huenemoerder et al. 


Robrade & Schmitt 


2007: 


Argiroffi et al. 



201 ll ). These densities, considering the 
typical accretion rates and surface fill- 
ing factors, are compatible with predic- 
tions of shock-heated material, and are 
significantly higher than that of typical 
quiescent coronal plasmas at tempera- 



Ness et al. 


2004: 


Testa et al. 


2004) 



10 



cm 



over iGiidel &: Telleschil ( 120071 ) observed a 
soft X-ray excess in CTTS with respect to 
WTTS, compatible with the scenario of a 
further plasma component at a few MK 
produced by accretion. However other 
results are discrepant with predictions: 
the observed Lx of the high-density cool- 
plasma component in CTTS is lower than 
that predicted from t he accretion rate by 



more than a factor 1 (lArgiroffi et al.ll2009 



Curran et ajj 120111 ). leaving the coronal 
component the major contributor to the 
X-ray emission in CTTS; furthermore in 
the cool plasma of CTTS the density in- 
creases for increasing temperature, at odds 
with predicti ons based on a single accre- 
tion stream ( Brickhouse et al. 2010l ). Be- 
cause of these apparent discrepancies dif- 
ferent scenarios were proposed, suggest- 
ing that the high-density cool plasma in 
CTTS could be coronal plasma, confined 
into magnetic loops, that is so mehow mod- 



2007 ; Brickhouse et al.ll2010l : iDupree et al. 



20121 ). 

In addition to containing plasma at a 
few MK, the shock region is known to be 
associated with material at T ~ 10 4 K or 
more, significantly hotter than the sur- 
rounding unperturbed photosphere, as 
a consequence of the energy locally de- 
posited by the accretion process. This 
photospheric hot spot produces excess 
emission in the UV and optical band, 
which is often rotationally modulated be- 
cause of the very small filling factor of 
the accretion-shock region and because 
accretion streams are usually not sym- 



(Bouvier et al. 


1993 


Petrov et al. 


2001 


)• 



Herbst et al.l 11994 : 



ified by the accretion process (IGiidel et al. 



Therefore, if the ob- 
served high-density X-ray emitting plasma 
also originates in the accretion shock, then, 
its X-ray emission might display rotational 
modulation. Specifically plasma heated in 
the accretion-shock, observed in the X- 
rays, could display periodic variations in 
density, emission measure, average tem- 
peratures, absorption, and source optical 
depth, as a consequence of stellar rotation. 
First hints of accretion driven X-rays that 
vary becaus e of the stellar rotati on were 
provided by lArgiroffi et al.l (120111 ) for the 
star V2129 Oph. 

Understanding the origin of this high- 
density plasma is important, both for 
constraining the total amount of X-rays 
emitted in CTTS, and setting the en- 
ergy balance of th e accretion-shock region 
(ISacco et al.l |2010| ) . Eventually, a defini- 
tive confirmation that this plasma com- 
ponent is material heated in the accretion 
shock, would make its X-ray radiation an 
insightful tool to probe the physical prop- 
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erties (i.e. density and velocity) of the ac- 
cretion stream, and to measure the chemi- 
cal composition of the inner disk material 
({Drake et alJl200fSh . 

To search for such X-ray modulation ef- 
fects we planned and carried out X-ray 
monitoring of V4046 Sgr, a close binary 
CTTS system in which both components 
are actively accreting from a circumbinary 
disk (see § E]) • 

In this work we describe the first re- 
sults from an XMM-Newton Large Pro- 
gram (LP) focused on V4046 Sgr, based on 
time-resolved high-resolution X-ray spec- 
troscopy on timescales down to 1/10 of the 
system orbital period. To constrain the 
large-scale magnetic field and the accre- 
tion geometry, we also carried out a coordi- 
nated multi-wavelength campaign involv- 
ing photometry, spectroscopy, and spec- 
tropolarimetry of V4046 Sgr. 

In §[2] we summarize the project focused 
on V4046 Sgr, whose properties are de- 
scribed in § [31 Details of the data process- 
ing and analysis are reported in § HI The 
observing results are presented in § [5j and 
then discussed in § [6j 



2. The V4046 Sgr project 

The XMM-Newton observation of V4046 Sgr 
consists of a 360 ks exposure performed on 
2009 September 15-19 (Obs-id: 0604860201, 
0604860301, and 0604860401). This ob- 
servation is part of a quasi-simultaneous 
multi-wavelength campaign (optical pho- 
tometry with REM/ROSS, 2009 Septem- 
ber 1-30; optical spectroscopy with TNG/SARG, 
2009 September 10-17; optical spectropo- 
larimetry with CFHT/ESPADONS, 2009 
September 2-8), aimed at studying simul- 



taneously the properties of coronal plas- 
mas, stellar magnetic field structure, pho- 
tospheric spots (both cool spots and hot 
spots), and the accretion process. 

Here we present the results obtained 
with the XMM-Newton/RGS specifically 
aimed at searching for rotational modula- 
tion in the accretion- driven X-rays. The 
results of the entire observing campaign 
are presented in a series of papers describ- 
ing, among other results, the properties 
of the X-ray emitting plasma (A. Mag- 
gio et al., in preparation), maps of the 
large-scale magnetic field structure and 
accretion geometry as inferred from opti 



cal spectropolarimetry ( iDonati et al.ll2011 
G. A. J. Hussain et al., in preparation, 
S. G. Gregory et al., in preparation), vari- 
ations in the accretion process over a 
range of timescales (G. G. Sacco et al., 
in preparation), detection and identifica- 
tion of a distant com oving WTTS system 
dKastner etld~]l201ll ). 



3. V4046 Sgr properties 

V4046 Sgr is a close CTTS binary sys- 
tem composed of two solar-like mass stars 
(masses of 0.91 and O.SSM^radii of 1.12 



and 1.04 R IDonati et al.ll201ll ). separated 



by 8.8 R . The two components are syn- 
chronously rotating wit h a period of 2.42 d, 



in ci rcularized orbits ( iStempels fc Gahm 



200J). V4046 Sgr is estimated to lie at 



a distance of 73 pc ([Torres et al.l 120081 ) 



and it is viewed with an inclination of 35 c 



the line of sight, 


Stempels &; Gahm 


2004; 


Kastner et al. 


2008 


), with the orbital axis 



of the binary likely aligned with the indi 
vidua ! stellar rotation axes (IDonati et al. 
201lh . At an age of ~ 10 - 15Myr 
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([Torres et al.ll200a bonati et alfcpllf ) and 
classified CTTS, V4046 Sgr is still 

surrounded by a dusty, molecule-rich cir- 
cumbinary disk ( iRodriguez et al.l 120101 ) 
from which bot h the components are ac - 
tively accreting (IStempels fc Gahml 12004 ). 



Newton observation^. In this w ork we 



A previous Chandra observation ( iGiinther et 



20061 ) showed that V4046 Sgr has a cool 
plasma component (T « 2 — 4 MK) at high 
density (n e « 0.3 — 1 x 10 12 cm~ 3 ), inter- 
preted as material heated in the accretion 
shock. 

At the time of the XMM-Newton ob- 
servation, the spectroscopic optical mon- 
itoring demonstrated that both compo- 
nents were accreting with a constant rate 
of 5 x 10~ 10 M yr _1 (infer red from the 



analy sis of the Ca II IRT, iDonati et al. 
201 ll ). Both components displayed com- 
plex magnetic fi elds (average surfac e inten- 
sity of ~ 200 G, IDonati et alll201lh . signif- 
icantly weaker than that of young er solar- 
like CTTS (e.g. IDonati et alJboiOl ). These 
magnetic fields are not strong enough to 
disrupt local disks farther than 1 R* above 
stellar surface, thus the formation of cir- 
cumstellar stellar disks around each com- 
ponent, distinct from the circumbinary 



disk, may be possible (jde Val-Borro et al. 
201 ll ). The accretion process, based on 



Ca II IRT, did not show significant ro- 
tational modulation, suggesting that post 
shock material contributing to these lines 
is symmetrically distributed with respect 
to stellar poles. 

The optical monitoring campaign con- 
firmed the orbital/rotational period (2.42 d), 
and determined the conjunction and quadra- 
ture epochs at the time of the XMM- 



adopt the phase reference defined in lStempels &: Gahm 
( 12004J . HJD = 2446998.335 + 2.4213459 E, 
with phase 0.0 indicating the quadrature 
with primary receding). However our op- 
tical monitoring revealed a phase shift of 
alP-069 with respect to that ephemeris, with 
quadratures occurring at phases 0.93 and 
0.43, and conj unctions a t phases 0.18 and 
0.68 (bonati et aDl201lh . 



4. Observations 

The XMM-Newton observation of V4046 Sgr, 
composed of three observing segments of 
~ 120 ks each separated by gaps of ~ 50 ks, 
covered 2.2 system rotations. X-ray emit- 
ting material heated in the accretion shock 
is expected to have temperatures of a few 
MK at most. Therefore to search for X- 
ray variability possibly produced in the 
accretion shock we analyzed the XMM- 
Newton/RGS spectra, that contain emis- 
sion lines that specially probe the coolest 
plasma components. 

The RGS spectrograph, composed of 
two nominally identical gratings (RGS1 
and RGS2), covers the ~ 2 — 38 A wave- 
length range. The first order spectrum, 
embracing the 4 — 38 A band, has a res- 
olution FWHM of 0.06 A, while the sec- 
ond order provides a resolution FWHM of 
0.03 A in the ~ 2 - 19 A range. We ex- 
tracted RGS spectra using the standard 
RGSPROC task. Data were filtered discard- 
ing time segments affected by high back- 
ground count rates. The final net expo- 
sures of the three observing segments were 
of 115, 122, 120 ks, respectively. We then 



■"■The quadrature with primary receding occurred at 
2455078.199 HJD 
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applied the rgscombine task to add the 
RGS1 and RGS2 spectra of the same or- 
der. Totally 34800 and 8200 net counts 
were registered in the first and second or- 
der RGS spectra, respectively. 

We analyzed the RGS spect ra using the 
IDL p ackage PINTofALE v2.0 flKashvap fe Drak e ; 
2Q00J and the XSPEC v!2.5 flArnaud 
19961 ) software. We measured individual 
line fluxes by fitting simultaneously first 
and seconcH order RGS spectra. Fit pro- 
cedure was performed in small wavelength 
intervals (AA < 1.0 A). The adopted best 
fit function takes into account the RGS 
line spread function (determined by the 
matrix response fuction), and the contin- 
uum contribution (determined by adding a 
constant to the line emission, and leaving 
this constant free parameter in the 

fit). 



5.1. Time resolved RGS spectra 



To investigate variability on short timescales, 
we analyzed RGS spectra gathered in time 
intervals of ~ 25 ks (i.e. bins of 0.12 in 
rotational phase). Totally nine lines have 
fluxes detected at la level in all the time 



5. Results 

The RGS spectra collected during the 
entire observation (see details in A. Mag- 
gio et al., in preparation) indicate that 
the main properties of the X-ray emit- 
ting plasma of V4046 Sgr are similar to 
those observed during the previous C han- 
dra observation (IGunther et al.l 120061 ): the 
plasma at T ~ 1 — 4MK has high den- 
sity, n e ~ 10 11 — 10 12 cm -3 , as determined 
by the f/i line ratio of He-like triplets of 
N VI, O VII, and Ne 1X0. 



2 Second order spectrum was used only for lines con- 
tained in its wavelength range. 

3 The measurements of the Ne IX triplet was per- 
formed by including in the fit the Fe XIX line at 
13.52 A, that is anyhow weaker than the Ne IX 
lines. 



intervals. These lines, and their fluxes at 
different time intervals, are reported in Ta- 
ble [TJ Significant variability on the ex- 
plored timescales is observed for all the 
listed lines. 

To check for variations in the coolest 
plasma components we considered lines 
with peak formation temperature T max < 
5 MK among the lines reported in Table [TJ 
This sample of lines, named cool lines, is 
composed of: the Ne IX triplet (13.45, 
13.55, and 13.70 A), O VIII Ly/fl and 
Lya (16.00 and 18.98 A), O VII resonance 
line (21.60 A), and N VII Lya (24.78 A). 
Among the lines reported in Table [TJ the 
Ne X and Fe XVII lines stay out of the cool 
line sample, because of their T max higher 
than 5 MK. Therefore their flux likely in- 
cludes significant contributions from hot 
plasma. These two lines compose the hot 
line sample. 

To maximize the S/N of the coolest 
plasma emission we added the measured 
fluxes of the cool lines for each time inter- 
val. This total line flux, plotted in Fig. [TJ 
is variable and the observed modulation is 
clearly linked to the stellar rotation: the 
flux is higher near phases 0.0 and 0.5, i.e. 
quadrature phases, and lower near phases 
0.25 and 0.75, i.e. conjunction phases. To 
confirm this variability pattern we fitted 



4 This is blended with an Fe XVIII, that is however 
negligible because of the EMD and abundances of 
the X-ray emitting plasma. 
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Table 1 

Observed line fluxes at different time intervals 



rot. 
phase 



Ne X 
12.13 A 
6.3 MK 



Ne IX 
13.45 A 
4.0 MK 



Ne IX 
13.55 A 
3.5 MK 



Ne IX 
13.70 A 
4.0 MK 



Fe XVII 
15.02 A 
5.6 MK 



O VIII 
16.01 A 
3.2 MK 



O VIII 

18.97A 
3.2 MK 



O VII 
21.60 A 
2.0MK 



N VII 

24.78 A 
2.0 MK 



-^intl/scgl 
-Pint2/scgl 
-Pint3/segl 
^int4/ scgl 



F 



nt5/ s 
-Fintl/s- 



-^int2/seg2 
-^int3 / seg2 
-^int4/seg2 
-Pint5/scg2 
^intl/seg3 
-Pint2/seg3 
-Pint3/scg3 
-Fint4/seg3 
-^int5/scg3 



0.67 
0.79 
0.91 
0.03 
0.15 
0.50 
0.62 
0.73 
0.85 
0.97 
0.32 
0.44 
0.56 
0.68 
0.79 



34.0±5.2 
33.0±5.0 
31.7±4.7 
45.9±5.4 
40.9±6.6 
40.2±5.2 
54.3±5.5 
45.6±5.5 
42.5±5.2 
42.6±5.7 
61.6±6.0 
48.3±5.6 
43.3±5.3 
23.7±4.3 
39.1±5.8 



58.3±6.6 
76.0±7.0 
56.1±6.9 
65.7±6.9 
59.8±8.3 
71.0±7.4 
37.7±6.1 
43.1±6.3 
55.4±6.6 
75.1±8.2 
46.3±6.0 
68.5±7.2 
66.6±7.1 
36.9±5.9 
41.9±6.4 



21.1±5.2 
17.7±5.1 
34.3±6.2 
31.6±5.8 
22.0±6.8 
33.8±6.1 
30.3±5.7 
45.4±6.2 
31.7±6.0 
60.8±7.5 
26.6±5.4 
41.3±6.2 
20.4±5.3 
9.9±4.3 
21.4±6.1 



21.3±4.8 
16.8±4.3 
25.3±5.5 
24.9±4.9 
27.1±6.3 
37.3±5.6 
35.9±5.3 
20.2±4.9 
29.5±5.3 
39.7±6.1 
20.7±4.6 
32.5±5.1 
35.9±5.1 
30.7±5.1 
20.7±4.9 



14.9±3.2 
7.4±3.0 
14.0±3.2 
10.5±3.3 
11.6±3.6 
14.5±3.4 
16.9±3.6 
18.5±3.7 
17.4±3.4 
5.2±2.8 
5.8±3.4 
14.7±3.7 
15.1±3.4 
4.7±3.0 
13.8±3.5 



15.3±3.3 
12.8±3.7 
17.6±3.7 

9.8±3.3 
22.3±5.0 
10.9±3.7 
13.1±3.7 
20.4±4.0 
10.6±3.4 
18.8±3.7 

9.7±3.7 
12.8±3.7 
16.6±3.9 
11.9±3.1 
13.7±3.7 



90.6±8. 

81.8±8. 
101.0±9. 
118.2±9. 

52.4±8. 

95.6±8. 

92.0±8. 
101.9±8. 
112.6±9. 
121.3±9. 
135.8±9. 

98.3±8. 
112.0±9. 

69.6±7. 

75.0±8. 



9 18.9±8.9 
1 17.0±9.3 
1 42.5±11.7 

7 48.3±11.0 
9 30.6±12.7 

4 33.2±11.8 

5 31.2±10.1 
4 19.0±7.2 
4 42.6±9.6 

6 24.1±9.7 

8 35.0±9.1 

7 53.3±12.3 
1 39.9±10.4 
7 24.9±9.9 

9 26.0±11.9 



67.6±9.4 

69.5±9.4 

89.9±11.0 

62.3±9.2 

52.3±10.8 

119.1±10.7 
64.3±9.3 
72.5±9.5 
77.0±10.4 

107.4±11.6 
69.4±10.0 
76.7±10.0 
94.6±10.9 
43.4±7.8 
56.4±10.0 



Note. — For each line, the Table head reports: ion, wavelength, and maximum formation temperature. F inti / Scg j 
refers to the line fluxes measured in the i— th interval of the j— th observing segment. For each time interval, the listed 
phase corresponds to the central time of the bin. Line fluxes are in 10 -6 phs -1 cm -2 . Errors correspond to la. 



Table 2 
Best fit parameters 



Set name 


Line set 


P a (d) 


A b (%) 




(Xred,l) d 


(Xred,2) C 


cool lines 


Ne IX+O VIII+O VII+N VII 


1.22 + 0.01 


23 + 2 


1.36 + 0.05 


2.91 


8.93 


cool line subset 


VIII+O VII+N VII 


1.21 + 0.02 


22 + 3 


1.35 + 0.06 


2.31 


5.50 


cool line subset 


VII+N VII 


1.22 + 0.03 


30 + 5 


1.30 + 0.05 


0.65 


2.89 


cool line subset 


Ne IX 


1.22 + 0.02 


25 + 4 


1.40 + 0.06 


2.33 


5.12 


hot lines 


Ne X+Fe XVII 


1.12 + 0.03 


14 + 5 


1.61 + 0.10 


3.63 


3.46 



a Period. 
b Amplitude. 

c First epoch of maximum flux after observation start, since HJD 2455089. 
d Rcduccd x 2 obtained with a sinusoid as best-fit function (4 free parameters). 
e Reduced % 2 obtained with a constant as best-fit function (1 free parameters). 
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Fig. 1. — Total flux of the cool line set versus time. The set of cool lines is composed 
of: Ne IX triplet (13.45, 13.55, and 13.70 A), VIII Lya and Ly/3 (16.00 and 18.98 A), 
VII resonance line (21.60 A), and N VII Lya (24.78 A)). Horizontal error bars represent 
the time-bin width. Dotted line marks the best-fit sinusoidal function. Orbital/rotational 
ph ases are computed a ccordi ng to the ephemeris HJD = 2446998.335 + 2.4213459 E defined 
in IStempels fc Gahml (120041 ) . Vertical dashed lines (dark gray) indicate quadrature and 



conjunction epochs, with the corresponding schematic views of the system plotted above 
(white and gray circles represent the primary and secondary components, respectively). Time 
intervals adopted for extracting spectra corresponding to low and high phases are marked by 
the vertical bands (light blue and light red for the low and high phase, respectively). 
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Fig. 2. — Total flux of the hot line set versus time. The set of hot lines is composed of: Ne X 
Lya line at 12.13 A and Fe XVII line at 15.02 A. Dotted line marks a sinusoidal function with 
the same period, phase, and relative amplitude obtained from the best fit of the total flux 
of the cool line set. The hot lines do not show rotational modulation, unlike the cool lines, 
see Fig. [TJ suggesting that their variability is associated with coronal plasma variability. 
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these observed flux variations with a sinu- 
soid plus a constant. We left all the best- 
fit function parameters (period, phase, am- 
plitude, and the additive constant) free to 
vary. We obtained a best-fit period of 
1.22 ± 0.01 d, and an amplitude of 23 ±2% 
with respect to the mean value (Table [2]). 
The inferred period is exactly half the ro- 
tational period of the system. As guessed 
maximum and minimum phases occur ap- 
proximately at quadrature and conjunc- 
tion, respectively. To check whether this 
observed modulation is effectively linked 
to the cool plasma emission, and not to a 
given line emission, we performed the same 
fit by separately considering the total flux 
obtained from different and independent 
cool line subsets. In all the inspected cases 
(see Table [2]) we found the same periodic 
variability (period, phase, amplitude). 

We checked whether this modulation 
is present also in the emission of hotter 
plasma by applying the same fit procedure 
to the total flux of the hot lines, Ne X and 
Fe XVII. Fit results are reported in Ta- 
ble [2J in this case the periodic modulation 
is not detected. The observed variability is 
instead likely dominated by hot (coronal) 
plasma. Figure [2] shows a comparison be- 
tween Ne X+Fe XVII line variability with 
modulation observed for the cool lines. 
The detected X-ray rotational modulation 
is also not visible in the EPIC lightcurves 
(A. Maggio et al., in preparation), even 
considering only a soft band. The substan- 
tial continuum contribution mostly due to 
the highly variable hot plasma likely masks 
the rotationally modulated signal. Hence 
we conclude that the observed X-ray line 
flux modulation is due to the high-density, 
cool plasma component. 



To understand the nature of the ob- 
served variability we searched for varia- 
tions in the average temperature by con- 
sidering ratios of lines originating from 
the same element. All the inspected ra- 
tios display significant variability, but are 
not correlated among themselves, and are 
not related to the rotational phase. We 
also searched for variations in the plasma 
density, probed by the f/i ratio of the 
Ne IX triplet. This line ratio is approxi- 
mately constant (f/i w 1, indicating n e ~ 
10 12 cm -3 ) during the entire observation, 
except for a lower value measured during 
the third interval of the second segment 
(f/i = 0.45 ± 0.13, corresponding to n = 
(5.2^fl) x 10 12 cm~ 3 ), and a higher value 
observed during the fourth interval of the 
third segment (f/i = 3l 2 '^, corresponding 
ton < 4x 10 11 cm -3 ). These variations ap- 
pear to be associated with episodic events, 
like clumpy accretion flows, and not with 
a rotational modulation effect. 

5.2. RGS spectra at different phases 

The total flux of the cool lines from 
V4046 Sgr displayed variations in time 
linked to the stellar rotation. To inves- 
tigate the differences of the X-ray emitting 
plasma corresponding to epochs of low and 
high fluxes of the cool lines, we added RGS 
data collected at the same phases with re- 
spect to the X-ray rotational modulation. 
We extracted two RGS spectra obtained 
by adding all the events registered dur- 
ing time intervals centered on maximum 
and minimum times, with duration of one 
fourth of the observed X-ray period (inte- 
gration time intervals are shown in Fig. [TJ. 
The two resulting low and high spectra, 
whose exposure times are 84 and 94 ks re- 
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Table 3 

Observed line fluxes at different phases 



Ion 


A 


T 


Ft 

± Low 


"High 




fAl 


(MK) 


(10 — 6 nlis~ 1 ™^ 2 1 




Ne X 


10.24 


6.3 


5.9±1.3 


4.1±1.2 


Ne IX 


11.55 


4.0 


4.1±1.6 


11.9±2.1 


Ne X 


12.13 


6.3 


46.2±3.0 


50.6±3.0 


Fe XXI 


12.28 


10.0 


4.3±1.6 


3.7±1.8 


Ne IX 


13.45 


4.0 


72.2±4.0 


75.3±3.9 


Fe XIX 


13.52 


10.0 


16.3±3.3 


9.4±3.3 


Ne IX 


13.55 


3.5 


27.3±3.1 


49.4±3.5 


Ne IX 


13.70 


4.0 


29.9±2.9 


41.7±3.1 


Fe XVIII 


14.20 


7.9 


5.1±1.6 


3.8±1.5 


Fe XVII 


15.02 


5.6 


13.6±1.8 


14.3±1.9 


O VIII 


15.18 


3.2 


4.3±1.7 


7.9±1.8 


O VIII 


16.01 


3.2 


19.5±2.1 


16.2±2.0 


Fe XVII 


16.78 


5.0 


11.4±2.1 


10.3±1.9 


Fe XVII 


17.05 


5.0 


12.7±2.4 


14.1±2.4 


Fe XVII 


17.10 


5.0 


5.8±2.6 


8.9±2.6 


O VIII 


18.97 


3.2 


90.5±4.7 


113.9±4.8 


N VII 


20.91 


2.2 


11.0±3.1 


7.0±2.9 


O VII 


21.60 


2.0 


24.8±4.0 


46.0±6.1 


O VII 


21.81 


2.0 


17.3±5.9 


24.4±6.3 


O VII 


22.10 


2.0 


12.2±3.9 


10.3±3.3 


N VII 


24.78 


2.0 


72.0±5.4 


106.6±6.0 


N VI 


28.79 


1.6 


25.0±4.5 


17.9±4.5 


N VI 


29.08 


1.3 


12.7±4.1 


7.5±3.9 


C VI 


33.73 


1.6 


25.3±5.2 


35.6±5.7 



Note. — Line fluxes measured during the high and low phases are listed. 
Flux errors correspond to la. 



Table 4 

Line flux ratios at different phases 





Line flux ratio 


Rlow 


Rm g h 


Ne IX 
VIII La 
N VII La 


(13.45 A) / Ne IX (11.55 A) 
(18.97 A) / VIII Ly/3 (16.01 A) 
(24.78 A) / N VII Ly/3 (20.91 A) 


6-5±?:i 


151J 1 



Note. — Ratio errors correspond to 68% confidence level. 
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Fig. 3. — RGS spectra corresponding to minimum (low) and maximum (high) phases, corre- 
sponding to exposure times of 84 and 94 ks respectively. For clarity reasons the two spectra 
are slightly smoothed, and the spectrum of the high-flux phase is shifted toward longer 
wavelengths by 0.1 A. 



spectively, are shown in Fig. [3j while the 
measured line fluxes, detected at la level 
in the two spectra, are listed in Table El 

We searched for differences in the low 
and high spectra to investigate how the 
emitting plasma properties vary between 
these two phases. The two spectra display 
significantly different photon flux ratios of 
N VII, O VIII, Ne IX lines, as reported in 
Table HI In principle, these line ratios may 
vary due to changes of absorption, plasma 
temperature, or source optical depth. In 
Fig. H] we plot the measured line ratios to- 
gether with the values predicted in the op- 
tically thin regime for different tempera- 
tures and different hydrogen column den- 
sities, Nr. 

Absorption can change line ratios be- 
cause, on average, lines at longer wave- 
lengths suffer larger attenuation for in- 
creasing Nr. The two N VII lines con- 
sidered here are an exception, because the 
absorption cross section of the interstel- 



lar medium has the oxygen K -shell edge 
(23.3 A, e.g. IWilms et al.ll2000l ) located be- 
tween their wavelengths, making the longer 
wavelength line, the Lya (24.78 A), slightly 
less absorbed than the Ly/3 (20.91 A). The 
two lines however suffer very similar ab- 
sorption, making the absorption effect of 
little relevance in the case of the N VII 
Lya/Ly/3 ratio (as can be seen from the 
upper panel of fig. IU where the curves pre- 
dicted for different Ah are very similar). 
Therefore any change in this line ratio, as 
that observed, is hardly explained in terms 
of Ah variability. Instead, an Ah decrease 
from the low to the high state might ex- 
plain the variation of the O VIII line ra- 
tios, but an opposite Ah variation should 
be invoked to justify the Ne IX variability 
(middle and lower panels of fig. HJ). All this 
findings indicate that the hydrogen col- 
umn density toward the source appears to 
be unchanged, and that line ratio variabil- 
ity is produced by a different mechanism. 
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This conclusion is supported by the similar 
fluxes between low and high spectra mea- 
sured for the two lines at long wavelengths 
(N VI at 28.8 A and C VI at 33.7 A), the 
most affected by absorption, and it is also 
confirmed by the full fledged analysis of the 
EPIC data presented in A. Maggio et al., 
in preparation, where Nr is found to vary 
by only a factor 2 over the whole observa- 
tion around a mean value of 3 x 10 20 cm -2 
(i.e. logiV H = 20.5). 

The three explored ratios depend also 
on temperature, because of the different 
energy of the upper levels of the two elec- 
tronic transitions considered in each ratio. 
In this respect, the three ratios do not vary 
consistently. In fact, a temperature de- 
crease from the low state to the high state 
could explain the increasing N VII and 
O VIII line ratios, but not the variation 
Ne IX lines. The derivation of the plasma 
model (see A. Maggio et al., in prepara- 
tion) is beyond the scope of this work, but 
we anticipate here that the EMD does not 
appear to vary enough between the two 
phases to justify the observed variations 
of the line ratios. Moreover the average 
plasma temperature (logT ~ 6.6 — 6.7), 
together with the measured A%, indicates 
that, in some phases, line ratios are not 
compatible with the optically thin limit, 
irrespective of the nature of their variabil- 
ity. 

Optical depth effects can change line 
ratios because each line optical depth 
is directly proportional to the oscillato r 
strength of the transition (jActonl Il978l ). 
Therefore, if optically thin emission does 
not apply, transitions with very different 
oscillator strength may suffer different at- 
tenuation/enhancement, with stronger ef- 



25 r 
20 1 
15 
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Fig. 4. — Lya to Ly(3 photon flux ratio 
for the N VII and O VIII H-like ions and 
13.45 A to 11.55 A for the Ne IX He- 
like ion, vs plasma temperature. Horizon- 
tal bands indicate values measured during 
high (red) and low (blue) phases. Black, 
dark gray, and light gray lines represent 
predicted values for different absorptions 
iVn, with labels reporting the correspond- 
ing logA^H value, and with thicker lines 
marking curve portions compatible with 
the observed ratios. 
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fects occurring in lines with higher oscilla- 
tor strengths. In the inspected ratios the 
lines with higher oscillator strength are 
the Lya line of N VII and O VIII, and 



the 1 3.45 A line of Ne IX (e.g. iTesta et al. 
20071 ). Non negligible optical depth is 



expected in strong X-ray resonance lines 
produce d from shock-heated plasma in 
CTTS flArgiroffi et all l2009h . The ob- 
served variable ratios might indicate that 
some lines are affected by a changing op- 
tical depth. Since the expected attenu- 
ation/enhancement with respect to opti- 
cally thin emission depends on the source 
geometry and viewing angle, stellar rota- 
tion can produce periodic changes in the 
line opacity, and hence in the observed 
ratios. However, once again the three ra- 
tios do not vary in the same direction, with 
N VII and O VIII ratios being higher in the 
low state, whereas the ratio of the slightly 
hotter Ne IX lines results higher in the high 
phase. 

Summarizing we stress the significant 
variations observed in line ratios between 
low and high phases. The origin of these 
variations remains unclear. If these varia- 
tions were due to changes in plasma tem- 
perature or absorption, a coherent behav- 
ior would be expected for the three ratios, 
and this is not the case. Opacity effects in- 
stead can operate in a more complex way, 
provoking both line enhancements or re- 
ductions, depending on the source geom- 
etry and viewing angle. This hypothesis 
is therefore the most intriguing, especially 
considering that in some phases line ratios 
are discrepant from the value expected in 
the optically-thin limit. 



6. Discussion 

The main result of the time resolved 
spectral analysis of the X-ray emitting 
plasma from V4046 Sgr (§ [5]) is that 
the high-density plasma component at 
3 — 4MK is rotationally modulated with 
a period of half the system orbital pe- 
riod, with maximum and minimum phases 
occurring at quadrature and conjunction 
epochs, respectively. The observed X- 
ray rotational modulation indicates that 
this high-density plasma component is not 
symmetrically distributed with respect to 
the stellar rotational axes. We also found 
that strong emission lines from this plasma 
component provide some indications of non 
negligible optical depth effects, and that 
the periodic modulation appears to be as- 
sociated with variations in the source op- 
tical depth, as evidenced by the significant 
variations in line ratios sensitive to opti- 
cal depth observed between low and high 
phases. 

The strongest X-ray emission lines, 
produced by shock-heated material in 
CTTS, are expected to have non-negligible 
optical depth due to the high density 
and typical size of t he post-shock region 
(lArgiroffi et al.l 120091 ) . Moreover the opti- 
cal depth should vary if the viewing geom- 
etry of the post shock region changes. 

Hints of non-negligible optical depth 
observed in the strongest X-ray lines of 
V4046 Sgr indicate that the high-density 
plasma is mostly concentrated in a com- 
pact portion of the stellar surface, as pre- 
dicted for the post-shock material. More- 
over, the variability of the optical depth 
can be naturally explained with the chang- 
ing viewing geometry of the volume occu- 
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pied by the high-density plasma during 
stellar rotation. This scenario requires 
plasma confinement by the stellar mag- 
netic fields. 

We observed an X-ray period of half 
the system orbital period, as already 
observed for accretion indicators (e.g . 



Vrba et al. 1993; Kurosawa et al. 



2005) 



and X-ray emission (jFlaccomio et al.ll2005f ) 
from some CTTS. That could be ex- 
plained, in the case of V4046 Sgr, by differ- 
ent scenarios. If the X-ray emitting plasma 
is located only in one of the two system 
components, then a period of half the ro- 
tational period is observed when there are 
two accretion-shock regions on the stellar 
surface at opposite longitudes, or there is 
only one accretion-shock region and the 
maximum X-ray flux is observed when the 
base of the accretion stream is viewed side- 



ways ( lArgiroffi et al.ll201ll a configuration 
that occurs twice in one stellar rotation). 

Considering the system symmetry and 
the accre tion geometry pre v iously sug- 
gested by IStempels fc Gahml (120041 ) , it is 
conceivable that both components pos- 
sess similar amounts of high-density cool 
plasma. In this scenario the half period 
can be naturally explained assuming that 
the location on each stellar surface of this 
plasma, compact and not azimuthally sym- 
metric with respect to each stellar rotation 
axis, is symmetric for 180° rotations with 
respect to the binary rotation axis. 

The simultaneous optical monitoring 
campaign indicated that the two compo- 
nents have similar accretion rates, vali- 
dating the assumption that the two com- 
ponents possess similar amounts of high- 
density plasma. However, the optical ac- 
cretion spots, probed by Ca II IRT, did not 



show rotational modulation (IDonati et al. 
201ll ). Therefore accretion regions emit- 
ting Ca II should be symmetrically dis- 
tributed with respect to the stellar poles. 
This scenario, different from that obtained 
from the X-ray data, could be reconciled 
considering that X-rays are likely produced 
only by a fra ction of the entire accretion- 
shock region (ISacco et al.ll20101 ). 

In conclusion, our XMM-Newton/RGS 
data of the V4046 Sgr close binary system 
have shown for the first time the rotational 
modulation of X-ray lines characteristic 
of a cool, high-density plasma corotating 
with the stars. This strongly support the 
accretion-driven X-ray emission scenario, 
in which the high-density cool plasma of 
CTTS is material heated in the accretion 
shock. It moreover suggests that the ac- 
cretion flow is channeled by magnetic field 
lines anchored on the stars, along small 
magnetic tubes. This is consistent with 
the general framework of magnetic accre- 
tion, but brings new insights into the ac- 
cretion mechanism in close binary systems 
of CTTS. 

This work is based on observations ob- 
tained with XMM-Newton, an ESA sci- 
ence mission with instruments and contri- 
butions directly funded by ESA Member 
States and NASA. C.A., A.M., and F.D. 
acknowledge financial contribution from 
the agreement ASI-INAF 1/009/10/0. 
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